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It is shown that high-energy features of the cosmic ray spectrum, in particular the kink around 4 PeV and 
the corresponding change in spectral index, may be explained from interactions between highly energetic cosmic 
protons and relic Big Bang antineutrinos, if the latter have a rest mass of about 0.4 eV/c^. This explanation 
is supported by experimental data from extensive air-shower experiments, and in particular by the observation 
(Fly's Eye) of a second kink around 300 PeV, and by the abrupt change in the chemical composition of the cosmic 
ray spectrum that occurs at that energy. Both facts follow naturally from our theory, which predicts additional 
verifiable features of the cosmic ray spectrum in the few-PeV region, e.g. an abrupt decrease in the p/a ratio. 
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1. Introduction 

In the 66 years that have passed since Pauli 
postulated the existence of neutrinos, a great deal 
has been learned about these evasive particles. A 
large variety of experiments have been performed 
on neutrinos produced by nuclear reactors, by the 
Sun, by particle accelerators, in the Earth's at- 
mosphere and in a Supernova explosion. How- 
ever, the most important (and in many respects 
the most interesting) source of neutrinos remains 
unexplored until today: The Big Bang. In my 
contribution to this conference, I will talk about 
the relic Big Bang neutrinos, I will suggest a pro- 
cess through which they might reveal themselves, 
and I will confront you with some experimental 
facts that seem to be strongly in support of this 
process. 

2. Relic neutrinos 

According to the Big Bang model of the evolv- 
ing Universe, large numbers of (electron) neutri- 
nos and antineutrinos have been around since the 
beginning of time. In the very first second, when 
the temperature of the Universe exceeded 1 MeV, 
the density was so large that the (anti-)neutrinos 
were in thermal equilibrium with the other par- 
ticles that made up the primordial soup: pho- 
tons, electrons, positrons and nucleons. Photon- 
photon interactions created e+e^ pairs, which an- 
nihilated into photon pairs. Interactions between 



(anti-)neutrinos and nucleons turned protons into 
neutrons and vice versa. 

This leptonic era came to an end when the den- 
sity dropped to the point where the mean time be- 
tween subsequent v interactions started to exceed 
the lifetime of the Universe, ~ 1 second after the 
Big Bang. The production of neutrons through 
DeP reactions had already stopped slightly ear- 
lier, since it became energetically impossible. The 
(anti-)neutrinos decoupled from the other parti- 
cles in the Universe and the n/p ratio determined 
the primordial He/H ratio. 

Since that moment, the wavelengths of the 
(anti-)neutrinos have been expanding in pro- 
portion to the size of the Universe. Their 
present spectrum is believed to be a momentum- 
redshifted relativistic Fermi-Dirac distribution: 

"<"> = SSSTklTT <^> 

characterized by a temperature T^, — 1.95 K, 
since {T^/T^f = 4/11 and T^ ^ 2.726 ± 0.005 K 
Q . The present density of these Big Bang relics 
is estimated at ~ 100 cm~'^, for each neutrino fla- 
vor. That is nine orders of magnitude larger than 
the density of baryons in the Universe. 

In spite of this enormous density, relic neutri- 
nos have until now escaped direct detection. The 
single most important reason for that is their ex- 
tremely small kinetic energy, which makes it dif- 
ficult to find a process through which they might 
reveal themselves. 



3. How to detect relic neutrinos? 

Let us imagine a target made of relic i>eS and 
let us bombard this target with protons. Let us 
suppose that we can tune this imagined proton 
beam to arbitrarily high energies, orders of mag- 
nitude beyond the highest energies that can be 
reached in our laboratories. Then, at some point, 
the proton energy will exceed the value at which 
the center-of-mass energy of the p — Ve system 
exceeds the combined mass of a neutron and a 
positron. Beyond that energy, the inverse /3-decay 
reaction 



P + Ve 



(2) 



is energetically possible. 

The threshold proton energy for this process 
depends on the mass of the D^ target particles. If 
this mass is large compared to the 10~^ eV kinetic 
energy typically carried by the target particles, 
this may be considered a stationary-target prob- 
lem, and the center-of-mass energy of the p — 9^ 
system can be written as 



EpTTi^ > 1.7 • 10^^ 



(eV)^ 



This process will thus take place when 



EpieY) > 



15 



1.7-10 
TOi/(eV) 



(4) 



(5) 



In our Gedanken experiment, this threshold 
would reveal itself through a decrease in the frac- 
tion of beam protons that traversed the target 
without noticing its presence, as Ep is increased 
beyond the threshold. Is there any experimental 
evidence for a phenomenon of this type? 

4. The cosmic knee 

Obviously, the target needs to be huge in order 
to see any effect of the described process. The 



largest target we can imagine is the entire Uni- 
verse, which is indeed filled with VgS. This tar- 
get is continuously traversed by protons, and by 
other atomic nuclei which together constitute the 
cosmic rays. 



„ 10 

> 5 

ai 

"V 2 
u 

i ' 

ftq 0.5 

-a 



■s 



0.2 



0.1 



"1 1 



TTTTtii — 1 iiiiiiii — I 1 1 1 mil — I nimn i miiiiii 



. ^4U^ 




! iiiiii| I I iiTTn| I iMiiii| t: 



O Akeno 
■ Fly's Eye 
♦ Haverah Park 



ml I . I NlJ I Muml I limill ill I M I I lllJIltlil 




1012 



1014 1016 lolS 

E (eV/nucleus) 



1020 



since nrii, <C rup. When the experimental mass 
value of the proton (938.272 MeV) is substituted 
in Equation S and the requirement is made that 
Ecni > mn + me (940.077 MeV), this leads to 



Figure 1. The all-particle cosmic-ray spectrum 

Figure H shows the spectrum of these cosmic 
rays, compiled by the Particle Data Group |^. 
The spectrum falls extremely steeply with energy. 
In general, it is well described by a power law 

S - -•■ («) 

with n w 2.7 for energies below 1 PeV. In order 
to display characteristic features of this spectrum 
which would otherwise be hard to discern, the dif- 
ferential energy spectrum has been multiplied by 
£;2-7 jj^ Figured. The steepening that occurs be- 
tween 1 PeV and 10 PeV, where the spectral index 
n changes abruptly from 2.7 to 3.0, is known as 
the knee of the cosmic ray spectrum. Ever since 
its discovery, this feature has been the subject of 
intense interest. 

We notice that this knee exhibits exactly the 
features that we expected to see in our Gedanken 
experiment: The particle flux suddenly starts to 
decrease when the threshold is passed. Therefore, 
we postulate the following hypothesis: 

The change of the spectral index in the all- 
particle cosmic ray spectrum at an energy of ~ 
4 • 10^^ eV is caused by the onset of the reaction 
p + i>e — > n + e'^, which becomes energetically pos- 
sible at this point. 

This hypothesis necessarily implies (Equation 
1^) that the mass of the electron neutrino equals 



0.4 eV. Also, the knee would have to be an exclu- 
sive feature of the proton component of the cosmic 
ray spectrum, if the hypothesis were correct. Be- 
yond 4 PeV, one would thus expect to see a grad- 
ual drop in, for example, the p/a event ratio, and 
this could be experimentally verified. 

5. A second knee?! 

If protons interact with the relic background 
neutrinos, other cosmic ray particles may as well. 
The equivalent reactions in which a particles are 
dissociated in collisions with relic neutrinos and 
antineutrinos 



a + Ve -^ 3p + n + e 
a + De —^ p -I- 3n -I- e^ 



(7) 
(8) 



have Q-values of 27.5 MeV and 30.1 MeV, re- 
spectively. The threshold energies for these re- 
actions are larger than the threshold energy for 
reaction (||) by factors of 61 and 66, respectively. 
The cosmic ray spectrum in the range around the 
"knee" was measured in detail by the Akeno ex- 
tensive air-shower experiment ^J . From a fit to 
their experimental data, the kink was located at 
log Ep (eV) = 15.67. 

If the above explanation of the origin of the 
knee is correct, a kink in the spectrum of cosmic 
a particles should thus be expected at an energy 
that is larger by a factor of 61-66: log Ea ~ 17.5 
{Ea « 3 ■ 10^^ eV). It is remarkable that sev- 
eral extensive air-shower experiments that have 
studied the cosmic ray spectrum at these high 
energies have reported a kink in the area around 
log-E = 17.5. The Fly's Eye detector, which 
obtained the largest event statistics, observed a 
change in the spectral index from 3.01 ± 0.06 for 
energies < 10^^'^ eV to 3.27 ± 0.02 for energies 
1017.6 ^ E < lOiS-5 eVg. The Haverah Park de- 
tector also reported a kink at log_E — 17.6, with 
the spectral index changing from 3.01 ± 0.02 to 
3.24±0.07§. 

The differential energy spectrum measured 
with the Fly's Eye detector is shown in Figure 
|p. The contents of the energy bins have been 
multiphed by E^ in this figure (compared to E"^-"^ 
in Figure p. The abrupt change of the spectral 
index n from 3.01 to 3.27 becomes very clear in 
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Figure 2. Experimental results obtained with the 
Fly's Eye detectorM]. See text for details. 



this way. 

The Fly's Eye detector is also capable of mea- 
suring the altitude at which the shower maximum 
is located. Figure 0a shows the amount of mate- 
rial (air) penetrated by the shower particles when 
they reach the shower maximum, Xmax, as a func- 
tion of the shower energy. It is remarkable that 
this distribution exhibits an unmistakable kink 
at approximately the same energy as the energy 
spectrum itself. This indicates a clear change in 
the chemical composition of the cosmic rays at 
that point. 

The dashed curve in Figure Ep represents the 
results of a simulation modeling showers initiated 
by protons [Q. The equivalent curve for showers 
induced by a particles can be derived from this 
in a straightforward manner. Two effects have to 
be taken into account. Both effects cause X^ax 
to be smaller for as than for protons of the same 
energy. 

The first effect concerns the interaction length 
(Aint) which for as is shorter than for protons. 
The interaction length of protons in air amounts 
to 90 g/cm^. Since the total cross section for 
interactions of highly energetic ions with atomic 



number A scales with A^^^, Aint scales with 
^-2/3^ This gives Aint = 55 g/cm^ for as in air. 

The second effect concerns the particle multi- 
plicity. In the first nuclear collision, the energy 
carried by the incoming a-particles is distributed 
to a larger number of secondary particles. In the 
subsequent shower development, the shower max- 
imum is thus reached earlier than in case of show- 
ers induced by protons of the same energy. The 
effect of this can, in first approximation, be esti- 
mated by considering the ion-induced shower as 
a collection of A separate showers, each with an 
energy Ea/A. From the slope of the dashed line 
in Figure 0a, one may estimate a-induced show- 
ers in air to reach their maximum earlier than 
proton- induced ones, by 30 g/cm^, counting from 
the first nuclear interaction. 

These (somewhat oversimplified) considera- 
tions thus lead to Xmax values for a-induced reac- 
tions in the atmosphere that are 65 g/cm^ smaller 
than the values for protons of the same energies, 
as indicated by the dotted curve in Figure go. 

The experimental data shown in Figures 0a and 
g& are in complete and detailed agreement with 
the expectations following from our hypothesis. 
These expectations can be summarized as follows: 

• At energies below 4 • 10^^ eV, the chemical 
composition of the cosmic rays is constant. 
Protons are the most dominant component. 
The average shower maximum follows the 
dashed line (Figure ^), as required by the 
well-known logarithmic energy dependence. 
The spectral index n — 2.7. 

• At 4 ■ 10^^ eV, the threshold for reaction 
(^ is crossed. The spectral index changes 
from 2.7 to 3.0. Beyond the threshold, 
protons gradually disappear from the spec- 
trum, which becomes more and more dom- 
inated by as. As E increases, the shower 
maximum thus moves away from the dashed 
line to the dotted line in the Xmax plot. 

• For energies around 10^^ eV, as have re- 
placed protons as the most dominant com- 
ponent. The spectral index is still 3.0 

• At 3 • 10^^ eV, a new threshold is crossed. 



this time for reactions (8) and (9). The 
spectral index changes from 3.0 to 3.25. Be- 
yond the threshold, as gradually disappear 
from the spectrum. In that process, pro- 
tons carrying a significant fraction of the 
energy are produced. This leads to a grad- 
ual enrichment in protons. Therefore, in 
the Xmax plot, the position of the shower 
maximum moves back from the dotted line 
to the dashed line, as the energy increases. 

At 10^^'^ cV, a new threshold is crossed, 
known as the ankle of the cosmic-ray spec- 
trum. At this energy, the gyromagnetic ra- 
dius of a proton starts to exceed the galac- 
tic radius and therefore extragalactic pro- 
tons may reach the Earth. The spectral 
index changes back to its original value of 
2.7. The chemical composition continues to 
change in the same sense as before, i.e. , a 
further gradual enrichment in protons. 
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Figure 3. The expected average depth of the shower 
maximum in the atmosphere (a) and the expected 
spectral index (fe), as a function of the energy of the 
cosmic rays. Experimental data were obtained with 
the Fly's Eye detector U. See text for details. 

This sequence of effects is graphically shown in 
Figure O, which is an extension of Figure and 
includes the knee area. The experimental data 
are in excellent agreement with these predictions. 



6. How is this all possible? 

We now turn our attention to the million-dollar 
question: How can the process that is the basis of 
our hypothesis play such a significant role, given 
its extremely small cross section? 

The cross section for the inverse process, 9^ 
scattering off protons, was measured to be pro- 
portional to the available energy in the center-of- 
mass system, i. e. the energy above the Q-value of 
1.805 MeVil: 



a = 0.34(£;c,„-Q)10 



-«cm2 



(9) 



with the energies expressed in units of MeV. 
Therefore, the cross section rises almost propor- 
tional to the energy of the incoming protons in 
the energy range between 10^® eV (cr « 10^"'"'^ 
cm^) and 10^^ eV {a « 10~^^ cm^). For a target 
density of ~ 100 cm~^, these cross sections trans- 
late in mean free paths of IQ^'^-^^ cm, or average 
lifetimes of iQ^^^^i years, i.e. 9-11 orders of 
magnitude longer than the age of the Universe. 
If this is all there is, the high-energy cosmic ion 
spectra could thus never have been affected at a 
significant level by the hypothesized process. 

However, it is important to realize that, with a 
mass of 0.4 eV, the relic PgS would be nonrelativis- 
tic {kT ~ lO^"' eV). Typical velocities would be 
of the order of 10^ m/s in that case[0, less than 
the escape velocity from the surface of the Sun. 
Such neutrinos may be expected to have accumu- 
lated in gravitational potential wells. WeilerS 
recently estimated that the density of relic neu- 
trinos in our own galaxy would increase by four 
orders of magnitude (compared to the universal 
density of 100 cm~^) if their mass was 1 eV. 

Locally, this effect could be much more spec- 
tacular. Extremely dense objects, such as neu- 
tron stars, black holes or active galactic nuclei 
(AGN) could accumulate very large numbers of 
relic neutrinos and antineutrinos in their gravi- 
tational field. Let us consider, as an example, a 
typical neutron star, with a mass (M) of 3 • 10'^*' 
kg and a radius of 10 km. Even at a distance (r) 
of one million kilometers from this object, the es- 
cape velocity is still considerably larger than the 
typical velocity of these relic neutrinos: 700 km/s. 



The concentration of relic neutrinos in such a 
local potential well is ultimately limited by the 
Pauli principle, which limits their phase-space 
density to Agi,h~^ , where g^, denotes the num- 
ber of helicity states and h Planck's constant |Q. 
Since the escape velocity scales with r~^/^, the 
maximum neutrino density is proportional to 
r^'^/^, and reaches values of O(10^^/cm^) near 
the surface of this neutron star. If the source of 
the potential well has a larger mass, the achiev- 
able neutrino densities increase proportional to 
M^/"^. In the "neutrino atmosphere" surround- 
ing an AGN with a mass of 1O^M0, the lifetime 
of a 10^^ eV proton could thus even be reduced 
to only a few years. 

The example of a neutron star was chosen since 
objects of this type may well play a role in acceler- 
ating protons to the very high energies considered 
here. Neutron stars usually rotate very fast and 
exhibit very high magnetic fields (typically ~ 10* 
T). When the magnetic axis does not correspond 
to the rotation axis, the changing magnetic fields 
in the space surrounding the neutron star may 
give rise to substantial electric fields, in which 
charged particles may be accelerated to high en- 
ergies 0. The strong magnetic fields, in turn, will 
tend to confine the accelerated particles to the 
region surrounding the neutron star. When the 
charged particles are electrons, the acquired en- 
ergy is emitted in the form of synchrotron radia- 
tion, a characteristic signature of many pulsars. 

However, for protons this type of energy loss 
only plays a significant role at energies that are 
13 orders of magnitude ((Tp/Te)**) larger than for 
comparable electrons. For example, a 10^^ eV 
proton describing a circular orbit with a radius of 
1,000 km loses per orbit about 0.6% of its energy 
in the form of synchrotron radiation. To put this 
in perspective, that is about a factor of five less 
than the fractional energy loss experienced by 100 
GeV electrons in CERN's LEP ring. 

It is, therefore, not inconceivable that a sta- 
tionary situation might arise in which protons 
and heavier ions, initially attracted by gravita- 
tional forces, and accelerated to high energies in 
the electric fields mentioned above, end up in 
a closed orbit around the dense object (neutron 
star, black hole, AGN) in which energy gained 



and synchrotron losses balance each other. In 
that case, the ions could spend a long time in the 
dense neutrino atmosphere surrounding the grav- 
itational source and collisions with relic neutrinos 
might become a realistic possibility. 

The idea that an AGN might act as the engine 
that produces cosmic rays is not at all new. If 
this engine works as a synchrotron, as described 
above, there are two ways in which the trapped 
accelerated ions can escape to other regions of 
the galaxy: Through collisions with other nuclei 
or through collisions with relic neutrinos. In the 
first case, the reaction products acquire trans- 
verse momentum and are scattered out of the ac- 
celeration plane. However, if a proton collides 
with a relic f'e, the resulting neutron has approx- 
imately the same four-vector as the original pro- 
ton. The escaping neutron thus remains in the 
acceleration plane. That also holds for the pro- 
ton into which this neutron decays after a while 
(100 years or so!). If the Earth was not located 
in the AGN's acceleration plane, these protons 
would thus never reach us, in contrast with the 
protons scattered in nuclear collisions. The mea- 
sured proton spectrum would thus exhibit a kink 
at the threshold for the pDe reaction. 

7. Conclusions 

The high-energy cosmic ray spectrum exhibits 
some intriguing features that can be explained 
in a coherent manner from interactions between 
cosmic protons or a particles and relic DeS if the 
latter have a restmass of 0.4 eV/c^: 

• Kinks at 4 ■ lO^^ gy and 3 • 10" eV, which 
would correspond to the thresholds for the 
pPe and ai^e reactions. 

• The energy separation between these kinks 
is consistent with the difference between the 
Q- values of the pDf. and aDe reactions. 

• The change of the chemical composition at 
3-10^^ eV, where as start to disappear from 
the spectrum and the protons reappear. 

If collisions with relic neutrinos were indeed re- 
sponsible for these features, a large concentra- 
tion of such neutrinos should be present in the 



vicinity of the source of the high-energy cosmic 
baryons. The required conditions could be met if 
these baryons were accelerated in the electromag- 
netic fields surrounding an AGN. 

There are various ways to test this hypoth- 
esis experimentally. First, one could explicitly 
measure the chemical composition of the cosmic- 
ray spectrum in the energy region from 10^^ to 
10"'^^ eV. Second, one could measure the shower- 
maximum distributions for energies below the 
Fly's Eye lower limit to see if the predicted kink 
near 4 PeV (Figure Ha) corresponds to reality. 

It should be emphasized that the predictions 
derived from our hypothesis are completely at 
variance with those derived from the "Standard 
Model of Cosmic Ray Physics" 0. In this 
model, the cosmic rays are the result of parti- 
cle acceleration in the shock waves produced in 
supernova explosions. This model predicts, for 
example, a maximum energy proportional to the 
nuclear charge Z of the particles. In this model, 
the "knee" would thus be a feature of the high-Z 
component of the cosmic rays, while we predict it 
to be a feature of protons and protons alone. 

Measurements of the high-energy cosmic pro- 
ton and a spectra would thus not only constitute 
a crucial test of our hypothesis, but also check 
some key predictions of this "Standard Model" . 
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